Coulomb interactions of massless Dirac fermions in graphene; pair-distribution 
functions and exchange-driven spin-polarized phases. 
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The quasi-2D electrons in graphene behave as massless fermions obeying a Dirac- Weyl equation 
in the low-energy regime near the two Fermi points. The stability of spin-polarized phases (SPP) 
in graphene is considered. The exchange energy is evaluated from the analytic pair-distribution 
functions, and the correlation energies are estimated via a closely similar four-component 2D electron 
fluid which has been investigated previously. SPPs appear for sufficiently high doping, when the 
exchange energy alone is considered. However, the inclusion of correlations is found to suppress the 
spin-phase transition in ideal graphene. 

PACS numbers: PACS Numbers: 73.43.-f,73.50.-h,73.22.-f 
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Introduction- The carbon atoms in graphene form 
a quasi-two-dimensional (Q2D) honeycomb lattice and 
contribute one electron per carbon to form an unusual 2D 
electron system (2DES) with a massless Dirac-fermion 
dispersion near the Fermi pointsQ, Q|. Graphene and 
related materials (e.g, nanotubes, fullerenes) have be- 
come a mine of novel technologies and a new paradigm 
for studying various aspects of physics 3] , including cos- 
mological models on honeycomb branes, superconductiv- 
ity on bi-partite lattices^ and nanotubes 0, Hubbard 
models^, spin-phase transitions^, and other aspects of 
strongly correlated electrons Q. 

The hexagonal Brillouin zone has two inequivalent 
points K=(l/3, 1/^/3) and K'=(-l/3, 1/^3), in units 
of 2n/ao, where ao is the lattice constant. The simplest 
tight-binding model with nearest-neighbour hopping t is 
sufficient to describe the low-energy regime where the 
valence and conduction bands (w and it*) have linear 
dispersion near the K, K' points, with zero bandgap. 
The graphene 2D electron system (G2DES) is nominally 
"half-filled" , with the it* band unoccupied, and has spin 
and valley degeneracies, with a Berry phase associated 
with the valley index 0]. 

The vanishing of the density of states and the effective 
mass of the 2D electrons near the Fermi points suggest 
that the Coulomb interactions of the massless fermions 
remain strong, unlike in the usual Fermi-liquid picture. 
This also means that perturbation approaches have to 
be treated with great caution. The Coulomb interac- 
tions may induce a gap between the it and it* bands, or 
a lifting of the sublattice (valley) degeneracies, or sta- 
bilize spin-polarized phases (SPP) in preference to the 
unpolarized state. Such SPPs in GaAs/AlAs 2DES, pre- 
dicted to appear at low coupling (^2 — 4) when pertur- 
bation methods are used, get pushed to high coupling if 
non-perturbative approaches were used(~ 26), as dis- 
cussed in Refs. HUGH EH- The two- valley 2D electron 
system (2v-2DES) does not show a SPP when full non- 



local non-perturbative calculations are used, presumably 
because these methods properly account for the direct 
interactions which are three times as many as the ex- 
change interactions^^- The exchange and correlation 
energy E xc in the 2v-2DES of the Si/Si02 system was 
calculated from the distribution functions in Ref . ^3] , us- 
ing the classical-map hyper-netted-chain (CHNC) tech- 
nique, accurately recovering the Quantum Monte-Carlo 
(QMC) results even in the strong coupling regimes |l3^. 
CHNC provides the pair-distribution functions (PDFs) 
gij(r) as a function of the coupling strength. Then 
E xc is evaluated via a coupling constant integration. 
The method has been successfully applied to a vari- 
ety of problems including the 2DES 14], the 2v-2DES 
in Si-MOSFET devices^, and the thick quasi-2DES in 
HIGFET structures 0. However, a direct calculation of 
exchange and correlation in graphene involves and 8x8 
matrix of two-component PDFs due to the spin and val- 
ley indices as well as the presence of tt, it* bands. Hence 
in this study we consider the exchange energy E x via an 
analytic evaluation of the non-interacting PDFs, and in- 
directly estimate the correlation energy E c using results 
for the spin-polarized four-component 2v-2DES, exploit- 
ing the analogies between the two systems. 

The non-interacting PDFs of G2DES, viz., 0?.(r) in- 
volve two components, the first being a Bessel function as 
in the ordinary 2DES, while a second, associated with the 
cosine of the angle of e-e scattering, involves Bessel and 
Struve functions, as we show below. We find that there 
are stable SPPs in an exchange-only approach. How- 
ever, including the correlation energy using the 2v-2DES 
data stabilizes the disorder-free G2DES in the unpolar- 
ized state. 

Model for the graphene 2DES -The kinetic energy near 
the K points is given by a Dirac- Weyl Hamiltonian of the 
form: 



H k = v F (p x T z a x + PyO-y) 
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FIG. 1: (Color online)Linear dispersion bands near a K point 
where the n* and 7r bands cross. In (a) we show a doped 
unpolarized system with equal occupation of the up-spin (a u ) 
and down-spin states. In (b) the polarized system has only 
electron carriers. In (c) both electron and hole carries occur. 
This is the only possibility if doping is zero (hp = and 



Here t z = ±1 defines the degenerate valleys, and a x , a y 
denote the x and y Pauli matrices that act in the space of 
the two atoms in each unit cell. The 7T, it* bands of spin 
and valley degenerate states (Fig.QJ show a linear disper- 
sion E = ±vphk. This form requires a cutoff momentum 
K c such that the number of states in the Brillouin zone 
is conserved. That is, if Aq is the area per carbon, then 
Kl = 47r(l/Ao). The electron density N c at half-filling 
is 1/Aq, with Ao = a.Qy/3/2, since one it electron of arbi- 
trary spin is provided by each carbon atom. The Fermi 
velocity vp — taoy/3/2 is thus the slope of the linear dis- 
persion, with vf ~ 5.5 eVA. If the G2DES is embedded 
in a medium with dielectric constant erj, then we define 
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(2) 



This is the ratio of a typical Coulomb energy to the hop- 
ping energy and hence is usually taken as the Coulomb 
coupling constant of the G2DES. This plays the same role 
as the r s parameter in electron-gas theory of nonrelativis- 
tic finite-mass fermions. The usual r s is not available for 
G2DES since the effective mass m* is zero and there is 
no effective Bohr radius. The coupling constant g° is 
maximized if eo is unity, and consistent with this case we 
assume g 
treat it as tunable 



- 2.75 for our G2DES studies, and do not 



The 4-component envelope-eigenfunctions of the ki- 
netic energy term are made up of two-component func- 
tions U = 0,e 4 * fc ), U' = (e^"» 6 ) and O = (0,0) where 
4>k is the angle of the vector k in the 2-D plane. Thus 



F^(r) = (2A)^(U 1 0) TX , 
FK(r) = {2A)^{0,U') TXa 



(3) 
(4) 



Here b = ±1 is a 7r*, tt band index, (• • -)t indicates the 
transpose, and \cr is the spin function. Then, using v — 
1,2 as a valley index, the Coulomb interaction may be 
written in the form: 



H, = 
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a k,Di ,bl ,CTl a p + q,D 2 M iCTl a P>«2,b3 ,C 2 a k+q,l>i ,bi ,<T 2 (^) 

Here a + , a are electron creation and annihilation opera- 
tors and V q = 2-ne 2 1 (t^q) is the 2D Coulomb interaction. 
The phase factors introduce a novel cos(6 l ) contribution 
where is the scattering angle, not found in the usual 
jcllium-2DES. The resulting form of the exchange energy 
per Carbon is: 



E x /E u 
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xkp 1 + b ^ 2 S2i^_ nfe (fc)n b 2|g (p) 
|k-p| 

Here we have introduced the intrinsic coupling constant 
g° and the energy unit E u = v F k c . Here k c = K c /^/2 — 
•y/(47rn c ) is based on the electron density per spin species, 
n c = N c /2 = 1/(2Aq). The above form of the exchange 
energy can be reduced to an evaluation of a few ellip- 
tic integrals!). The normal "half-filled" G2DES can be 
doped with electrons or holes; but it is easy to show that 
symmetry enables us to limit to one type of doping. How- 
ever, given a system with an areal density of N$ dopant 
electrons per valley, with n$ = N$/2 per spin, the carri- 
ers in the spin-polarized system could be electrons only, 
or both electrons and holes, as shown in Fig. ^ f° r the 
7r* and 7r bands at one K point. The intrinsic system 
with ns = can be an unpolarized state, as in Fig.^a), 
or spin-polarized state with electrons and holes, as in 
Fig. ^c) . Such exchange-driven systems have been stud- 
ied by Peres et al.||, while the correlation effects have 
not been considered. Since the correlation energy terms 
out-number the exchange terms by 3:1, it is imperative to 
include the correlation energies. The calculation of cor- 
relation energies is always more difficult than exchange 
energies. Peres et al. treated g° as an externally tun- 
able parameter for forcing SPPs. Here we evaluate the 
exchange energy E x at the intrinsic g°, from the non- 
interacting PDFs, and indirectly evaluate the correlation 
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FIG. 2: (Color online)The Bessel-like and Struve-like non- 
interacting, parallel-spin PCFs h b (r) and h s (r) for the un- 
polarized doped system. The bands are numbered as in 
Fig. 0b). The anti-parallel non-interacting PCFs are zero. 



energy E c from the four-component 2v-2DES with the 
same coupling strength (r s = g ) and spin-polarization. 

Pair- distribution functions of the G2DES- Although 
we are dealing with a four-component system (2- valleys, 
2-spin states), as seen from Fig.Q] we need to consider the 
redistribution of electrons and holes among the tt* and 
7T bands when comparing the energy of spin-polarized 
states with the unpolarized state. From Fig. we see 
that the e-e interactions at a given valley can be con- 
structed from (i) interactions with a tt*(7 u band of up- 
spin electrons of density n u , filled to k u , (ii) a n*<Jd set of 
electrons or a iro~d spin-down holes, of density rid, filled 
up to kd (iii)the tt<j u band, with electron density n c , filled 
to k c and (iv)the irud band, density n c , filled to k c There 
will also be similar inter- valley terms. Each term in this 
4x4 matrix, denoted by Qij{r) where i, j = 1 • • -4, will 
have two components associated with those in U and 
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(6,e 4 ^). Thus g i3 {r) 



WJ ~ 9ij(r),9ij(r), where the 
superfixes "&, s" indicate that the noninteracting forms 
are Bessel-function like, and Struve-function like, respec- 
tively. The Struve form arises from the cos(#) terms in 
the Coulomb interaction. The numbering scheme of the 
matrix is shown in Fig. 0b). Thus, defining the pair- 
correlation functions (PCFs) Hij(r) = Qij{r) — 1, or its 
components hij{r) = gij{r) — 1, we have: 
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h °ij S ( r ) = 1 
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Here Jo, Ji are Bessel functions, while Hq and Hi are 
Struve ii-functions. Also, in [JqHi — Jii?o]i the func- 
tions are evaluated at the argument fc^r. The wavevec- 
tors ki = y/^Trrii) are for each component i, of density n,. 
We show typical noninteracting PCFs for a doped, un- 
polarized case as in, Fig. 0a), with the doping fraction 
ns/n c = 0.2. In CHNC, the exchange-hole is mapped 
exactly into a classical Coulomb fluid using the Lado 
procedure 0|- The figure shows that the exchange- hole is 
strongly reduced by the presence of the cos(#) term which 
has been averaged into the Struve-like PCFs h s (r). When 
the Coulomb interaction is included, the cos(#) term has 
a similar mitigating effect and exchange-correlation in 
the G2DES is actually considerably weaker than in the 
corresponding 2v-2DES which may be realized in a Si 
MOSFET. The CHNC calculation for the 2v-2DES for 
the conditions stipulated in Fig.|21show that the correla- 
tion energy is only about a third of the exchange energy 
at r s = g . This and other model calculations allow us 
to conclude that a good approximation is to ignore the 
E c contribution arising from the Struve-like distribution 
functions and estimate the E c only from the Besscl-likc 
functions which are the only forms appearing in the 2v- 
2DES. This justifies our use of the 2v-2DES for the cor- 
relation energy, while the E x is exactly evaluated. 

The kinetic and exchange energies.- When the dop- 
ing per valley is N$ = 2ns, the total number of elec- 
trons per valley is N t = N c + Ng. Also, using the 
i = 1,2,3,4 notation of Fig. H£b), we se t n i = n ui 
n-2 = nd, n 3 = n 4 = n c . Hence the spin polariza- 
tion s = n u — bdUd, where the band index bd = — 1 
for holes. The degree of spin-polarization £ = s/N t . 
The composition fractions, inclusive of the valley index 
v =1,2 are x V i = m/2Nt- We note that fc_F = ^(2-Kns), 
k u = \/{2ir{ns + s)}, kd = ^/{2tt\hs - s\}. The exchange 
energy E x (n$, £) can be written as: 



E x (ng,0/N t = (N t /2) 



2nrdr 



(8) 



It is implied that the Struve-like component in Q v ,v,ij 
where v label the valleys, is summed with the appropri- 
ate bibj band ± factors. Only a sum over the compo- 
nents in one valley is needed in evaluating the exchange. 
The calculation of both exchange and correlation involves 
exactly the same formula, but with a further integra- 
tion over a coupling constant A included in the PDFs 
Qv,v',ij.\- The above formula for E x is made more ex- 
plicit below. Thus, the total kinetic and exchange energy 
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Ekx= K.E+£ , 2 ,, for the case (b) of Fig.^can be written 
in terms of np, n u , nd and n c as in Eq. [HJ or in terms of 
hp, k u , kd, k c and A as: 



Ekx (C) 
A 



(in 



Mkl + kj) 



(9) 



(27T) 



G?V/4)(7r/2)[fc*Wu(r)+fc3W a2 (r) 



2fc c 2 {fc^ 13 (r)+fc^ 14 (r)}] 

The kinetic and exchange energy, E tx (( = 0) of the un- 
polarized system shown in Fig. [Tfa), is obtained by set- 
ting kp = when k u = kd in Eq. |5| The energy dif- 
ference A.Efcx(C) is plotted in the upper panel of Fig. 13 
This equation reduces to Eq.(12)-(15) of Ref. |3| when 
the PDFS are replaced by their fc-space forms and ex- 
pressed in terms of elliptic integrals. We have done the 
calculations in r-space using PDFS, and in fc-space via 
elliptic integrals, as a numerical check. The evaluations 
using elliptic intergrals are numerically more stable. The 
results (Fig. EJ) show that stable spin-polarized phases 
appear for systems containing purely electron carriers, 
when the dopant fraction fd > 0.38. 

Effect of the correlation energy E c .- This requires the 
calculation of all the intra and inter-valley interacting 
PDFs of G2DES. Even using various symmetries, this 
involves calculating nearly two dozen PDFs as well as an 
integration over the coupling constant. We postpone this 
demanding task to future work. However, the E c (ns, C) of 
G2DES may be quite well approximated from the closely 
analogous four-component 2v-2DES. This correlation en- 
ergy/per electron, e(r s , £), where r s — <?°, is given in Eq. 
5 of Ref. [l^. I n transferring from the 2v-2DES to the 
G2DES we note that e 2 /e which is unity in 2v-2DES 
becomes g°vp in the G2DES. Also, the exchange term 
arising from the Bessel-like hn(r), and li22(r), calculated 
at the given coupling strength r s = g° and the spin po- 
larization £ arise from the same PDFs in both systems. 
This establishes the scale factor connecting the energy 
units of the 2v-2DES and the G2DES at any given spin- 
polarization. Given the validity of Eq. |H1 we adopt the 
same scale factor for converting the correlation energy 
E c of the 2v-2DES to the G2DES. The energy difference 
A-Efc xc (C)j inclusive of E c , between the polarized and un- 
polarized phases calculated using the above approach is 
shown in the lower panel of Fig. [3] The exchange-driven 
SPP seems to be suppressed to within the uncertainty 
of this calculation of AE^xc- Improved estimates would 
only strengthen this calculation further. 

In this work we have kept the Coulomb coupling 
fixed ato° = 2.75 typical of graphene, unlike in other 
studies[a 13 where the coupling strength g is taken as a 
tunable parameter, (in the spirit of Hubbard- model stud- 
ies). Even in one- valley 2DES, the SPP of low-order 
theories is pushed to g ~ 26 — 27. In the 2v-2DES, 
direct terms (usually known as the correlation terms) 
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FIG. 3: Upper panel- the energy difference AEtx, i.e., 
K.E+exchange, between the polarized and unpolarized 
phases, in units of E u — vpk c , as a function of the spin polar- 
ization £ and the dopent fraction fd~ns/n c . Electron-carrier 
systems, Fig.^b) are more stable than electron-hole systems, 
and show stable spin-polarized states. However, addition of 
the correlation energy (lower panel) makes the unpolarized 
state the most stable phase. 



predominate over exchange interactions, and the SPP is 
suppressed, as found in CHNC|I3 or QMC[l3|. Here 
too, we conclude that the inclusion of correlations sup- 
presses the exchange-driven spin-phase transition in the 
ideal graphene system. 
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